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Abstract —The downlink of a single cell base station 
(BS) equipped with large-scale multiple-input multiple-output 
(MIMO) system is investigated in this paper. As the number of 
antennas at the base station becomes large, the power consumed 
at the RF chains cannot be anymore neglected. So, a circuit 
power consumption model is introduced in this work. It involves 
that the maximal sum-rate is not obtained when activating all 
the available RF chains. Hence, the aim of this work is to 
find the optimal number of activated RF chains that maximizes 
the sum-rate. Computing the optimal number of activated RF 
chains must be accompanied by an adequate antenna selection 
strategy. First, we derive analytically the optimal number of 
RF chains to be activated so that the average sum-rate is 
maximized under received equal power. Then, we propose an 
efficient greedy algorithm to select the sub-optimal set of RF 
chains to be activated with regards to the system sum-rate. It 
allows finding the balance between the power consumed at the 
RF chains and the transmitted power. The performance of the 
proposed algorithm is compared with the optimal performance 
given by brute force search (BFS) antenna selection. Simulations 
allow to compare the performance given by greedy, optimal and 
random antenna selection algorithms. 

Index Terms-Massive MIMO, number of RF chains, antenna 
selection, circuit power consumption. 

I. Introduction 

A major challenge of designing future wireless networks is 
to support the huge increase in data traffic demand. Massive 
MIMO system is one important technology that will be 
adopted in the fifth generation (5G) of cellular networks m, 
O, Q), j4j. It is based on the installation of large number of 
antennas at the Base Station (BS). The growth of the number 
of antennas involves quasi-orthogonality between the users’ 
channels in consequence of the law of large numbers. Hence, 
a huge gain in spectral efficiency and energy efficiency can 
be obtained while using low complexity transmit and receive 
techniques such as Maximal Ratio Combiner (MRC), Max¬ 
imal Ratio Transmission (MRT) and Zero Forcing (ZF) 0, 

na, q. 

Few works have previously investigated the optimization 
of resource allocation in massive MIMO systems. In GO. 
the authors proposed an algorithm to compute the optimal 
power allocation among users in downlink massive MIMO 
systems with MRT precoder. However, a low complexity 
linear precoder is designed for downlink massive MIMO 
systems in 0. The beamforming vector is computed jointly 


with optimal power allocation in order to maximize the sum- 
rate. The proposed precoder archives high performance near 
Regularized Zero Forcing (RZF) and better than MRT with 
no matrix inversion. In iflOl . the authors investigate the uplink 
of multi-cell massive MIMO system with MRC receiver. 
An efficient algorithm is proposed to compute jointly the 
optimal training duration and power allocation. In HD, the 
number of RF chains is assumed to be limited. Hence, 
a joint antenna selection and user scheduling strategy is 
introduced for downlink massive MIMO systems. As the 
large number of antennas involves the growth of the Channel 
State Information (CSI), the authors in ( T2| assume limited 
CSI at the BS and efficient antenna selection algorithm is 
proposed. In CD. the authors propose a polynomial time 
algorithm to select antennas with maximum SNR. 

In massive MIMO systems, as the number of antennas 
is getting larger, the power consumed by the RF chains 
becomes not negligible D3, uni, m, D3. Consequently, 
the maximal sum-rate is not obtained when activating all 
the RF chains. Hence, the optimal number of RF chains to 
be activated that maximizes the sum-rate must be derived. 
In liT6ll . the authors investigate the optimization of the power 
allocation and the number of transmit antennas based on 
the asymptotic approximation of the capacity. They extend 
their work in na to multi-users and imperfect CSI case 
assuming MRT precoding. But, the optimization of the power 
allocated among users is still based on the average capacity 
over channel realizations. It will be more interesting if the 
optimization is based on the instantaneous throughput. As 
results, the optimal number of transmit antennas should be 
derived jointly with antenna selection strategy depending of 
the channel gain coefficients. In this work, we derive first 
analytically the optimal number of activated RF chains that 
maximizes the sum-rate averaged over channel realizations 
under the assumption of equal received power. Then, an 
efficient greedy algorithm is proposed to compute the close- 
to-optimal number of RF chains and to select the best subset 
of antennas that maximizes the sum-rate. The proposed 
algorithm allows to find the balance between the power 
consumed at the RF chains and the output transmitted power. 
Monte-Carlo simulations allow to compare the performance 
of the proposed algorithm with random and brute force search 
(BFS) optimal antenna selection. 


In this paper, all power variables are assumed to be unitless 
since they are normalized by the average noise power. [_-J 
denotes the floor function, [.] denotes the ceiling function, 
diag(p) is a diagonal matrix whose diagonal entries are 
the elements of the vector P, (x) + denotes max( 0 ,at), (.) H 
represents the Hermitian matrix, C^, is a binomial coefficient 
defined as Cfj = (jy-'s)!JV! an ^ Tr{.} denotes the trace of 
a square matrix. 

The rest of the paper is organized as follows. In Section II, 
the system model is presented. In Section III, the problem is 
formulated. In Section IV, the optimal number of RF chains 
that maximizes the average sum-rate is analytically calculated 
under the assumption of equal received power between users. 
Then, an efficient greedy algorithm is proposed to compute 
the number of RF chains that maximizes the sum-rate and to 
select the best antennas in Section V. The computational com¬ 
plexity of the proposed algorithm is calculated and compared 
to the optimal BFS algorithm in Section VI. Numerical and 
simulation results are shown and discussed in Section VII. 
Finally, we conclude and discuss the main results in Section 
VIII. 


II. System Model 

We consider the downlink of a single cell massive MIMO 
systems. The base station (BS) is equipped with a large 
number of antennas N serving I\ single-antenna users with 
N K. The channel gain is represented by complex matrix 
G = D 1 / 2 !!. The small scale fading is defined by complex 
matrix H = [hi, I12 ,Iik], where hk € C lxAr is the 
k th channel vector for user k, is assumed to be quasi-static 
Gaussian independent and identically distributed (i. i. d.) slow 
fading channel. Considering only path loss, the large scale 
fading is defined by the vector d = [d{ - a di a ...d~ K a ], where a 
is the path loss exponent and d k the distance between the BS 
and user k, D is a diagonal matrix defined as D = diag(d). 
The noise is assumed to be additive Gaussian white (AWGN) 
random variable with zero mean and unit variance. The vector 
P = [piP 2 ---Pk\ is assumed to be the distribution of power 
among users. We consider that BS knows perfectly the CSI. 
We consider the Zero Forcing (ZF) as beamforming strategy 
because it cancels the inter-user interference and it achieves 
high performance 0, (7). The beam forming matrix is 

expressed as 


W = 


H h {HH h )~ 1 

7 ? 


( 1 ) 


where 7/ the normalization factor defined as 77 = 

Hence, the received downlink vector of signals is expressed 
as 


y = 


£>5 

-dtag(^p) s + n, 

arj 


( 2 ) 


where s is the vector of transmitted data symbol and n is the 
vector of AWGN. a = y is a normalization factor. 


The signal received by user A: can be written as 


Uk = y/Pk J ^s k + n k . (3) 

arj 

The received SINR for user k is expressed as 


d, a pk 

SINR k = (4) 

a z rj z 

Hence, the sum-rate is expressed as 


K 

r = Y, 1o & 

k= 1 



d k >k \ 

cr 2 r) 2 ) 


(5) 


Since the number of antennas is large in the considered 
system, we assume a non-negligible circuit power consump¬ 
tion on, na. Let p c denotes the fixed power consumed 
at each activated RF chain (Digital to Analog Converter 
(DAC), mixer, frequency synthesizer, filter) and p max the 
maximal available power at the BS. Hence, the circuit power 
consumption constraint can be expressed as 


N 

Pout "F ^ \ O^n-Pc ^ Pm ax’) (6) 

n—1 

where a n is an antenna coefficient that is set to 1 if antenna 
n is activated and to 0 otherwise and p out is the output 
transmitted power given by 


K 

Pout = J2p k . (7) 

k =1 

It is to be noted that \p m ax/Pc\ represents the maximum 
number of RF chains that can be supported (assuming no 
transmission power) by the system due to the circuit power 
constraint. 


III. Problem Formulation 


The circuit power consumption involves that the maximal 
achieved sum-rate is not obtained when activating all RF 
chains. Hence, the optimal number of RF chains Y^n =1 °n 
that maximizes the sum-rate should be derived. Moreover, 
the aim is to find the optimal balance between the power 
consumed by the RF chains a n-Pc and the trans¬ 

mitted power p ou t. The number of RF chains should be 
optimized jointly with adequate antenna selection strategy. 
The transmitted power must be allocated optimally among 
users. Hence, the main problem can be formulated as 










K 


maximize 

R = 

p.cfn :n=l. .N 

K 

subject to 



k =1 


N 


E 

n =1 


r =J2 1o s2 1 


k= 1 


d~ k a Pk 


N 


a 2 r] 2 (ai,a 2 ,...,a n ) 


: Pmaxi 


a n G {0,1}, n = 1..N. 

( 8 ) 

It is to be highlighted that the term t) 2 in the objective 
function depends on the fading coefficients of the selected 
antennas. So, an antenna selection strategy is necessary to 
find the optimal number of RF chains that maximizes the 
sum-rate. Consequently, the problem becomes combinatorial 
with exponential complexity growth in N. It is known that, 
after selecting the best antenna, water filling is the optimal 
strategy for power allocation among users m- 


IV. Equal Received Power 

In this section, equal received power between users is 
assumed and the sum-rate averaged over channel realizations 
is approximated. Then, we derive analytically the optimal 
number of RF chains that maximizes the sum-rate averaged 
over channel realizations. The number of RF chains is given 
by 


N 

S = Y,«n. (9) 

n=l 

Note that since the sum-rate is averaged over the channel 
realizations, the antenna selection is meaningless. Since equal 
received power between all users is assumed, the transmitted 
power allocated for user k can be expressed as 


Pk 


Pmax S.p c <7“ 


K 


d~ k a 


( 10 ) 


The term rf is approximated in |5j, (6), J7J, i ffS) , when 

K, S —> oo 


rj 2 Tr{(iTiT ff ) -1 } K 

This approximation is validated by simulations in Section 
VII. Hence, the sum-rate averaged over the channel realiza¬ 
tions can be expressed as 


R = K. log 2 



0 


max 


S.p c ){S-K) \ 
K 2 ) ' 


( 12 ) 


In consequence, the optimal number of RF chains that 
maximizes the average sum-rate over channel realizations can 
be derived 


S* = 


where 


L0J if #(L</>J) > R (\H) or L'/’J = L Pmax/Pc\ 

|~ ( t>~\ otherwise 

(13) 


Pmax + A -Pc 

2p c 


< Pmax / Pc- 


(14) 


Finally, the power allocated for each user k is given by 


Pk 


Pmax - S*.p c (T 2 


K 


d k ° 


(15) 


However, the number of RF chains that optimizes the 
instantaneous sum-rate (for one channels realization) should 
be derived. It depends of the coefficients of the selected 
antennas. In consequence, the number of RF chains should 
be optimized jointly with adequate antenna selection strategy 
which will be the focus of the next section. 


V. Efficient Greedy Algorithm 

The aim of this work is to find the optimal number of 
RF chains and have an antenna selection that maximizes 
the sum-rate. The number of RF chains must be derived 
jointly with adequate antenna selection strategy. The antenna 
selection is a combinatorial problem with exponential growth 
in N. Moreover, the optimal antenna selection can be done 
with high complexity brute-force search (BFS) algorithm. 
So, we propose a low complexity greedy iterative algorithm 
that determines the number of RF chains maximizing the 
instantaneous sum-rate. At each iteration for a fixed number 
of RF chains, the best antenna n* that maximizes the sum- 
rate is determined among the set of non selected antennas A. 
Since, the cost function rj 2 is infinite for .S' < K as the ZF 
Precoder does not exist, we choose the first K antennas by 
optimal BFS. The best antenna should minimize the term rj 2 

argmax A f? = argmin A ?y 2 (ai, a 2 , ■■■, cuv) (16) 
Hence, the best antenna at each iteration can be derived as 


n* = argmin A 77 2 (ai,a 2 ,...,Q!Ar) (17) 

As the selected antennas are given, the optimal strategy for 
power allocation among users is given by Water filling lH9l 


- ( 1 _ °~y y 

Pk Un(2)M d~ a ) 


(18) 


where /j is the Water level. 

The convergence of the proposed algorithm is obtained 
when the instantaneous sum-rate starts decreasing (R$ < 
R S - 1 ). The convergence point exists since the sum-rate is 
a concave function, the second order derivative is negative 

j2 d 

< 0. Hence, the proposed algorithm allows to determine 
the number of RF chains, the selected antennas, the power 
allocated among users and the achieved sum-rate. 










Algorithm 1 Greedy Algorithm 
1: Initialization 

2: a n <— 0, n = 1 : N 
3: A <— {n, n = 1 : N} 

4: for S = K + 1 : \p m ax/Pc\ do 

5: Select the best antenna n* from A that minimizes r/‘ 2 

6 : OL n * i — 1 

7: A-^A\{n*} 

8: EfcL 1 (eMU - if) = Pmax - S.p c > find /X 

by bisection method 

9: Pfc , k = 1 : A', power alloca¬ 

tion 

10 : Break If Rs < Rs-i 

11: end for 


VI. Complexity Analysis 

In this section, the computational complexity order for 
both optimal brute force antenna selection and the proposed 
greedy algorithm are calculated. For the brute-force search 
optimal algorithm, the number of possible combinations of 
sets of the antennas is 


Hence, the optimal antenna selection can be obtained 
with very high complexity. However, the proposed efficient 
algorithm is polynomial time. The complexity order for both 
algorithms are presented in Table Q] 


TABLE I 

Computational Complexity of the proposed greedy and BFS 

ALGORITHMS. 



BFS 

Greedy 

Complexity order 

Cbfs 

Cgreedy 

N = 64 , K = 10 , lPmax/Pc] = 32 

0(8.10 16 ) 

0(6.10 5 ) 

N = 128 , K = 10 , lPmax/Pc] = 128 

0(3.10 3 °) 

0(4.10 7 ) 


VII. Numerical Results 

In this section, we present numerical results to show the 
performance of the proposed greedy algorithm compared to 
random antenna selection and optimal BFS algorithm. Monte- 
Carlo simulations are done to show the optimized number of 
RF chains and the maximal achieved sum-rate. 

Simulation parameters are summarized in Table HU The BS 
is equipped with 256 antennas serving 10 users. 


L Pmax I Pc J 

E C N-s+ 1 (19) 

S=1 

For the proposed greedy algorithm, the number of combi¬ 
nations of sets of the antennas is given by 

[Pmax /Pc\ 

E N-S+l (20) 

s=i 

For each selected set of antennas, the coefficient rj 2 is 
calculated. So, we need one matrix multiplication and one 
matrix inversion. The complexity order of the matrix multi¬ 
plication is 0{SK 2 ) and the complexity order for the matrix 
inversion is 0(K 3 ). 

Hence, the computational complexity order of the BFS 
algorithm is given by 


L Pmax /Pcj 

Cbfs = 0 | E C s n _ s+ 1 (SK 2 + K 3 ) 
s =1 


( 21 ) 


and the computational complexity order of the greedy algo¬ 
rithm is given by 


C, 


greedy 


= o 


\Pmax/Pc\ 

E 

S=1 


{N-S + l) (S K 2 



= 0{\jp max /p c \K 2 {K{N + l) 

+ 7} (N + 1 — K ) ( [Pmax/Pc\ +1) 

+ c(L P max /Pc\ + 1)(2| _p max / Pc J + i))) 
D 


( 22 ) 


TABLE II 

System Parameters. 


Symbol 

Description 

Value 

K 

number of users 

10 

N 

number of antennas 

256 

Pc 

fixed power RF chains 

0.05 

a 

path loss exponent 

3.7 


cell radius 

500 m 


Fig. Q] shows the maximal achieved sum-rate in function 
of the number of antennas N under the proposed greedy 
algorithm and optimal antenna selection. Due to the high 
complexity of BFS algorithm, simulation results for optimal 
antenna selection are presented for limited number of anten¬ 
nas N, three users K = 3 and p rnax = 1- It is clear that 
the achieved sum-rate increases as the number of antennas 
increases since it offers more diversity to the transmitter. The 
difference between the performance of the low complexity 
proposed greedy algorithm and the optimal antenna selection 
strategy is small. 

In Fig. [2] we see the achieved average sum-rate in function 
of the number of RF chains under the proposed greedy 
algorithm and random antenna selection for different values 
of total available power at the BS. First, the obtained results 
confirm the hypothesis that the maximal achieved sum-rate 
is not obtained when activating all RF chains. Then, the 
performance of the proposed greedy algorithm is significantly 
higher than the random antenna selection because the greedy 
algorithm allows to select the best antenna at each iteration. 

Fig. [3] shows the optimal number of RF chains to be 
activated that maximizes the average sum-rate in function 
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Fig. 1. Maximal achieved sum-rate under optimal and greedy antenna 
selection (K = 3 ,Pmax = !)• 


Fig. 3. Number of activated transmit RF chains under random and greedy 
antenna selection. 



Fig. 2. Average sum-rate under random and greedy antenna selection. 

of the maximal available power at the BS. We observe that 
the analytic expression of the optimal number of RF chains 
obtained under equal received power and random antenna 
selection fits with simulation results. Then, it is clear that 
the optimized number of RF chains obtained by the proposed 
greedy algorithm is less than the number of RF chains 
obtained under random antenna selection due to the adequate 
antenna selection strategy and the power is optimally dis¬ 
tributed among users. So, the proposed algorithm allows to 
use less number of transmitted RF chains and more power is 
available at the base station for transmitting. Hence, it allows 
to increase the achieved sum-rate. 

Now, we investigate the impact of the number of users K 
on the system performance. In Fig. [4] we see the maximal 
achieved sum-rate in function of the number of users under 
the proposed greedy algorithm and random antenna selection. 
So, the optimal number of users that maximizes the sum-rate 
can be observed for different values of that maximal available 
power Pm ax at the BS. The proposed greedy algorithm allows 
to transmit with less number of RF chains than random 


antenna selection. So, there is more available power for 
transmitting p ou t- Hence, more users can be served under the 
proposed greedy algorithm and the maximal achieved sum- 
rate increases significantly. 



Fig. 4. Impact of the number of user K on the system performance. 

Finally, Fig.0shows the optimized number of RF chains in 
function of the number of users K for different values of the 
maximal available power at the BS. The number of RF chains 
increases when increasing the number of users K since we 
need more RF chains for transmitting. Also, as the maximal 
available power p ma x increases, there is more flexibility to 
use more transmit RF chains. It is clear that the computed 
number of RF chains under the greedy algorithm is less 
than random antenna selection. In consequence, there is more 
available power for transmitting and the system performance 
is ameliorated. 

VIII. Conclusion 

The downlink massive MIMO system with non negligible 
circuit power consumption is investigated in this paper. 










































Fig. 5. Number of activated transmit RF chains in function of the number 
of users K. 


Since this work considers the circuit power consumption, 
the maximal sum-rate is not obtained when activating all 
the available RF chains. For this reason, we focus in this 
paper on optimizing the number of activated RF chains 
that maximizes the sum-rate. So, we derive analytically the 
optimal number of activated RF chains that maximizes the 
average sum-rate under equal received power. However, the 
number of RF chains that maximizes the sum-rate should 
be optimized jointly with adequate antenna selection. Hence, 
an efficient greedy iterative algorithm that optimizes the 
number of activated RF chains is proposed. It allows to 
find the balance between the power consumed at the RF 
chains and the transmitted power. Simulations are done to 
compare performance of the proposed algorithm to random 
and optimal antenna selection. Numerical results show the 
efficiency of the proposed algorithm. 
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